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Abstract 

We determine the propagation properties of a quantum particle in a d- 
dimensional lattice with hopping disorder, delta correlated in time. The system 
is delocalized: the averaged transition probability shows a diffusive behavior. 
Then, superimposed to the disorder, we consider a bias favouring the motion 
with a given orientation, as in the dynamics of flux lines in superconductors. The 
result is an effective Liouvillian for the density matrix, which is characterized 
by competition between single particle and pair hopping. In this case the tran- 
sition probability is determined in terms of excitonic motion, each exciton being 
extended along the bias direction. In the small bias regime the hopping disorder 
is almost uneffective along the Bragg lines of the Brillouin zone, where drift 
dominates. Elsewere the system undergoes diffusion. In the opposite regime we 
find the single-sided hopping spectrum, as expected from the bias term, but, due 
to the hopping disorder, this undergoes an abrupt change of sign at the Bragg 
lines. 



1 Introduction 

Various studies have been devoted to quantum propagation in disordered lattices, in- 
cluding site and hopping disorder. Here we study time dependent hopping: in general 
the adiabatic motion of a particle in a "hot" background. The motion of a charge in a 
rapidly fluctuating effective magnetic field belongs to this class: here we give a lattice 
version of the problem. Time-dependent fluctuations of the magnetic field have been 
considered by Aronov and Wolfle in studying the behavior of doped high-T c materials, 
close to the metal-insulator transition Q : their analysis was motivated by magnetore- 
sistence measurements pj in Bi 2:2:0:1 compounds. Tight binding hamiltonians were 
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also considered for the dynamics of flux lines in superconductors, a widely investi- 
gated topic, both at the experimental and at the theoretical level: see, e.g.,Q,||] and 
references therein. Columnar defects, artificially produced by energetic heavy ion 
radiation, have been used in order to pin flux lines and reduce dissipation. Greatly 
enhanced pinning has been obtained, e.g., in YBa2Cu^O'j crystals with aligned colum- 
nar defects, produced by Sri- ion radiation In the corresponding path integral 
description, the euclidean time is the the vortex line parameter and 

the horizontal coordinates of the columnar defects || define the lattice nodes. 
In a hollow cylindrical superconductor the longitudinal current creates a transverse 
magnetic field which forces the flux lines to tilt with respect to the vertical alignment. 
In the hamiltonian this translates into a term, linear in the momentum [jl0| |Tl| , and 
antihermitean. It can obviously be read as originating from an imaginary vector po- 
tential. This term explicitly breaks the space inversion symmetry: in fact the particle 
has different left and right hopping amplitudes (in a given direction). This is not to 
be confused with a chiral particle, characterized by single-sided but unitary propaga- 
tion. Motion with a preferred orientation arises in various non-quantum mechanical 
contexts, e.g. in population dynamics Jl2| , in the transport of passive scalars in fluids 
[0, in directed percolation. The nonhermitean hopping term has the effect of depin- 
ning the vortex lines, as shown by various authors: @,||,(l4),@. Here we first 
consider time-dependent hopping with no bias. Due to the averaging, the appropriate 
object to be studied, rather than the wave function, is the density matrix. In terms 
of it one reconstructs every transition probability. Our approach relies on a second 
quantization formalism, which proved to be very efficient in describing edge states in 
quantum hall systems Q . In the limit of fast fluctuations memory effects are canceled 
and the effective dynamics is described by a Liouvillian operator (see Section 2). We 
find that the quantum particle undergoes classical diffusion. In Section 3 we add the 
deterministic bias: the Liouvillian, in the second quantization formalism, takes then 
the form of the so-called pair hopping model hamiltonian |TiJ , . Pure diffusion is 
now always frustrated; we find excitonic states, propagating with a nontrivial dis- 
persion law. In Section 4 we summarize our results and compare them with related 
work. In Appendix 1 we derive a property for averages of time-ordered exponential 
operators; in Appendix 2 we show that the transition probability, in the hermitean 
case, can be obtained by resumming the ladder diagrams, i.e. that it coincides with 
the diffuson amplitude. 

2 Disordered lattice 

We start with a lattice hamiltonian with time-dependent hopping disorder, including 
an antihermitean term: 

H (t) = - [(u(x,n;t) +w(x,n;t))\x X x + 

X,fJ, 

+ (u*(x, n; t) — w*(x, fi; t)) \x + >< x\] 
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here \i — 1,2, ...d, d is the lattice dimension. We assume zero average gaussian 
coefficients; the various amplitudes are mutually independent, with correlators: 

< u(x, /i; t)u*{x', fjf; t') > = 5(t — t') ■ S(x — x') ■ 5^ ■ D v (x,fi\t) 
<w(x,fi;t)w*(x,fx;t)> = 5(t — t ) ■ S(x — x ) ■ • D w (x, fi;t) 

As previously announced, we deal with the density operator in the second quantization 
formalism. We introduce two, mutually commuting, fermi (or equivalently bose) 
operators a(x) and b(x), related with two independent copies of the system. Before 
averaging, they are associated with the retarded and advanced particle , and evolve 
independently with evolutors U and U* . We define the operator F, as the second- 
quantized evolution operator: its matrix elements in the 1+1 particle sector are: 

< 0\b(y) ■ a{x) ■ F(t,t') ■ a+(x') ■ b+(y')\0 >=< x,y\U{t,t') ® U*(t, t')\x', y' > (1) 

F is associated with the following two-particle hamiltonian: 

H(t) = -^[u{x^-t)-C{x^)+u*{ Xl ii-t)-C + {x^) (2) 

-I- w(x, fi; t) ■ B(x, /it) — w*(x, /it; t) ■ B + (x, /i)] 

C(x,n) = a + (x) ■ a(x + e M ) — b + (x + e M ) • b(x) 

B(x,n) = a + (x) ■ a(x + e M ) + b + (x + e M ) • b(x) 

Q{x) = a + {x)-a{x)-b+{x)-b{x) 

Using the fact that the disorder is ^-correlated in time, one can exactly perform the 
average of the Neumann series and reexponentiate the result, thus obtaining: 

<F(t,t')>=Texp[-sign(t-t') J dr ■ -ff e //(r)] (3) 

The effective Liouvillian H e f f has the form: 

H ef f{t) = 5 ^ *) ' M). C + (x, m)} - D w (x, w t) ■ {B{x, M ), B+(x, /x)}] , 

(4) 

here { , } denotes the anticommutator. The average has generated a quartic term, 
which couples the particle and the antiparticle. The nonlinearity can be easily handled 
in this case: in fact -ff e // can be written as a quantum spin hamiltonian Q, if one 
starts with the fcrmionic representation. One verifies that an angular momentum 
algebra is obtained from a(x) and 6(a;). The angular momentum is given by: 

J+(x) = a + (x)-b+(x) (5) 

2-J 3 (x) + l = a + {x)-a{x)+b + {x)-b(x) = N(x) 

If we define D± iX)fJ/ (t) — D u (x, /x; t) ± ■_D t0 (x, /j,; t); the Liouvillian turns into: 

H ef f(t) = -^^2[i-D +tX ^(t)-(j 1 (x)-J 1 (x + e tl ) + J 2 (x)-J2(x + e ll )) 

X,fJ. 

+ D- tX)li (t) ■ (4 • J 3 (x) ■ J 3 (x + e^) + Q(x) ■ Q(x + e M ) - l)] 
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The planar term, which describes pair hopping, is ferromagnetic. The vertical term, 
which counts the particles, turns from ferro to antiferromagnetic as the antihermitean 
disorder overcomes the hermitean one. The angular momentum operators commute 
with the charge operators Q(x). Similarly the total number N a of a-type particles, 
Ni, , and Q(x) commute with the hamiltonian. Obviously, as long as we are concerned 
with the density matrix, we are only involved in the N a = Nb = 1 sector, as made 
explicit in the matrix elements written in Eq.|l|. In a first class of eigenstates the 
particles are separated and do not propagate, since the hopping term acts only on 
doubly-occupied sites, i.e. on pure states of the form \x >< x\. Let us denote 
such eigenstates as localized. A second class is given by plane wave superpositions 
(magnons) of doubly-occupied sites (in the case of homogeneous disorder): 

J + (P)\0>= j^ m Y^j + (x)-exp( l p-x)\0> (|*v|<tt) (6) 
with eigenvalues : 

E(p; t) = 2 J2lD-At) ~ D+,»(t) ■ cos P^ ( ? ) 

The site transition probability is then decomposed in plane wave contributions: 

(| < x\U(t,t')\x' > | 2 ) = . . , [ dp- explip- (x — x') — sign(t — t') [ dr ■ E(p; r)l . 

[2n) (t J_ v J t , 

(8) 

The magnons are insensitive to any site potential: this rather unintuitive result 
depends on the delta-correlation of the hopping amplitudes as shown in Appendix 
A3. Hence, with time-independent correlators and hermitean disorder, the Liouvillian 
reduces to a lattice laplacian: it is no surprise then that diffusive (long range order) 
modes are the outcome of averaging over fast time- fluctuations. If the disorder has 
an antihermitean part, hopping rates locally break space inversion invariance, and 
make some diffusive modes unstable. In fact, the spectrum of the Liouvillian is no 
longer positive definite: from the minus sign in the exponent in Eq. ^, one has that 
the portion of the Brillouin zone inside the surface E(p) = becomes unstable. Let 
us comment on the connection with the problem of a d=2 particle in a magnetic 
field. In two dimensions, a fluctuating magnetic field orthogonal to the plane would 
be described by an hopping coefficient of the form: u(x, /i, t) — expii ■ 6(x, /i, t)) with 
9(x, fi, t) gaussian; in our model instead u is gaussian. In spite of this major difference, 
the present solution confirms a previous result on the motion in a fluctuating magnetic 
field, obtained in the continuum case |lq] . The effective motion of the quantum 
particle is in both cases classical diffusion. We are not able to understand wether 
this is a mere coincidence or is related with some general property shared by both 
approaches. 

3 Biased system with disorder 

We add now a deterministic asymmetric hopping, which describes a biased transport 
in a preferred direction p2l. The hamiltonian, for the particle- antiparticle system, 
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becomes: 

5[ 1 ^ia[exp{— k)a + (x + e M ) • a(x) + exp(k)a + (x) ■ a(x + e M )] 

^ ^ i7I Q;[exp(-A;)6 + (a; + e^) • b(x) + exp(k)b + (x) ■ b(x + e M )], (9) 

where k and a are real. Since in the perturbative series one can isolate the deter- 
ministic term and expand in the disorder term (see Appendix Al), the total effective 
Liouvillian L is simply: L — i ■ sign(t — t') ■ H e ff(t) + Hbias, where H e ff is given 
in Eq. ^. If we consider homogeneous hopping disorder, we recover a non-hermitean 
version of the so-called pair hopping model hamiltonian jl8| . Notice that here the pair 
hopping term is intrinsically dissipative. Neither the magnons nor the eigenstates of 
the bias term (free-particle states) are eigenstates. It is nonetheless possible to deter- 
mine two families of solutions, which can be regarded as the natural extension of the 
previously determined ones (localized and diffusive, respectively). The wave function 
f(x,y), in the two-particle sector (N a = Nb — 1), satisfies the eigenvalue equation 

a ■[e- k f(x-ej T ,y)+e k f(x + ej T ,y)-e- k f(x,y-ej T )-e k f(x,y + e JT )} (10) 

+ i ^2 D +(^) s x, y [fi x + e M>y + + f( x - e ^v - - 2D -(v).f( x ,y) = £■ f(x,y), 

where, with respect to the previous section's notation, we have £ = —i ■ E. Since the 
pair-hopping term vanishes on singly-occupied sites, two-particle states, separated 
in every other direction but the bias one, will only be acted by the single hopping 
term. This first class of solutions is the obvious extension of the formerly localized 
ones. In the second class there is no longer separation orthogonal to the bias. Let us 
proceed to the details of the solution. Upon writing / as a function of the baricentric 
coordinate R = (x + y)/2 and of the relative coordinate r = x — y, one easily identifies 
the eigenspace S°, which can be spanned in terms of the eigenfunctions: 

fs,n,B,o(R, r) = exp(iP ■ R + icfc- r^-) ■ [IL^f^ S R ^ R oJ, (11) 

where R a ,n are(d-l)-dimensional vectors, playing the role of degeneracy indexes. 
They are the projection of the baricentric and relative coordinates on the space Xj_, 
orthogonal to the bias. The eigenvalues then, up to a constant, coincide with the 
eigenvalues of Huas ■ 

£{Pp,Qii) = -4a-sin[(P/2) 1I -ik}-sin[cti T }-2iJ2 D -(v)- (12) 

A* 

One verifies that S° is eigenspace also upon adding a static disordered potential 
V(x). Spectrum and eigenfunctions will then reproduce the features discussed by 
Hatano and Nelson (l(| . To go over to the second class of solutions we first Fourier 
transform the equation: 

{£ - 2a ■ cos[(P/2 - qfc - ik] + 2a ■ cos[(P/2 + q) ¥ - ik])f(P, q) = (13) 
= -2^£>_(/x)/(P,q)+2^ j D + (/x)co S (P / ,)-^ f " dqf(P,q) 



H, 



bias 
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After integrating over q, we get: 



+7r dq' 



r 1 =(£/4)+i/2-J2D-(ii); z = a-sin[P T£ /2-ik] ;^^D+( M WP M ), 

where the integral is one-dimensional. The eigenvalues arc: 

£{P) = -2* • £ D_(aO + 4i£ • [1 - (^) 2 ] 1/2 (15) 

Notice that the solution is invariant under the symmetry P — ► — P , fc — > — k; out of 
the two branches of the square root, we take the one that, as a — > 0, goes into the 
spectrum of the unbiased case (Eq. |^). The wavefunction fp (P, q), associated with 
£(Po), has the form: 

fro(^) = , s%) + iMl _ (f]1/ r (16) 
In coordinate representation one finds an exponential behavior: / as ^ , where: 

C ± = ^-[[i-(f) 2 ] 1/2 ±i], CC+ • C- = -1) (17) 

and ± is to be chosen according with the condition |£±| < 1. 

The function / is divergent as r-p — > — oo, but when computing the transition 
probability between two sites one only needs r-p = both in / and in the solution of 
the transpose equation. We finally examine the case with £± lying on the unit circle. 

The condition \£±\ = 1 implies real, with absolute value smaller than one: 

this is verified only in the absence of asymmetry in the deterministic term (A: = 0), or 
when Pp = ±7r. We have an exciton, extended along /j, = JI , with global momentum 
Po and relative momentum p e f f : 

Peff = arm8 L.Jnm)ir/2] )ik = 0) ' (18) 

Peff = arcosi \— -) (Pjt = ±7r) 

y a ■ coshyky * 

The spectrum of the Liouvillian E(P) (Eq. [l5|) in d=2, with bias in the ~p = 1 direc- 
tion, and with hermitean, isotropic disorder ( D w (fi) — 0, D u (fi) — D), is exhibited in 
the Figures. The imaginary part Ej{P) of E(P) describes the reversible motion, with 
drift velocity v = VpEj(P). The real part Er(P) describes the irreversible motion: 
for k 7^ a region of instability at the center of the Brillouin zone appears. When a/D 
is small enough, i.e. when disorder dominates, drift is found along the Bragg lines 
P2 = ±Pi ± 7r (see Fig.l); out of such lines Pr(P) is practically zero and the exciton 
predominantly undergoes diffusion. Upon increasing a, the bias exponent k tilts the 
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plane Ej(P) in the Pi direction, thus enforcing the drift also at the center of the Bril- 
louin zone. This is exibited in Fig. 2 At very strong bias ( k large) one would expect the 
single-sided hopping to dominate. The situation is rather different (see Fig. 3). One 
easily verifies that E(P) w sign[(cos(pi)+cos(p2)]-{ceexp(k))-(cos(pi/2)+i-sin(pi/2), 
where the last factor is indeed the eigenvalue of the bias operator in the single-sided 
hopping limit. The effect of disorder is in the first factor: an abrupt change in sign 
of both the dissipative and reversible part of the spectrum ay the Bragg lines. Let us 
finally discuss the case k = 0, which describes a particle with deterministic anisotropic 
diffusion plus disorder. On qualitative grounds, everything goes as in Fig.l, but now 
Er{P) is always positive (no instability occurs). Reversible and irreversible motion 
are now completely separated, in dependence on P. Around the Bragg lines, we have 
only drift, apart from a constant damping factor; in the complementary region we 
find pure diffusion, with Ej{P) — (no drift). The drift region tends to broaden as 
a is increased. 



4 Conclusions 

We have discussed the motion of a particle over a lattice with fastly fluctuating hop- 
ping amplitudes: the model describes a massive quantum object coupled with a high 
temperature background. In our formalism the transition probability is written as a 
transition amplitude for a two-particle quantum system. This makes simpler the op- 
eration of averaging the probability over classical fluctuations. Before averaging, the 
two particles evolve independently, respectively forward and backward in time; after 
averaging, they interact and their motion becomes irreversible. Their effective hamil- 
tonian, which is simply the Liouvillian of the density matrix for the original system, 
has a quartic interaction, which can be readily put in the form of an Heisenberg hamil- 
tonian. We determined the steady states of the Liouvillian: with them one computes 
the density matrix and any single-particle transition probability (see Eqs.1,8). A first 
class of steady states has the two particles physically separated and localized. The site 
transition probability depends on double occupancy states: it is a sum of plane waves, 
evolving in time with a diffusive law. The generic behavior of such quantum systems 
is then diffusion, and this holds true even in the presence of a disordered potential, 
as shown in Appendix 3: since the hopping amplitudes are delta-correlated in time, 
they destroy the phase coherence of the wave function. Quantum interference effects, 
essential for localization, are thus absent. Diffusion in a quantum- mechanical system 
was found in the Harper model at its critical point |l^],(20). The present result, ob- 
tained from a lattice model, confirms a previous one, on the motion of a particle in 
a rapidly fluctuating magnetic field, derived in the continuum case Jig] ,. In Section 
3 we added a deterministic, anisotropic bias, enforcing a favoured orientation along 
a given direction; as already illustrated, this term arises quite naturally in describing- 
tilted vortex motion in superconductors. The Liouvillian takes then the form of the 
so-called pair-hopping model hamiltonian (i.e. it includes both single-particle and pair 
hopping, and the two terms do not commute). In our context the coupling constants 
are complex, since we are mixing reversible and irreversible motion. Two classes of 
steady states can be found. In the first class the two particles do not interact, pro- 
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vided that their wave packets do not overlap in the plane orthogonal to the bias. On 
such states the hopping disorder has practically no effect, and the interesting physics 
is the depinning transition, as described by Hatano and Nelson. In the second class 
the two particles form an exciton, extended along the bias direction. The dispersion 
law is a nontrivial function of the exciton momentum. For small enough bias, the 
Brillouin zone splits into a diffusion-dominated and a drift-dominated part, the latter 
lying around the Bragg lines. The site transition probability, which adds over the the 
exciton contributions, is the sum of two parts: essentially reversible evolution around 
the Bragg lines, and irreversible diffusion with no drift in the complementary region. 
Notice that this is different from a mere sum of the two types of motion, since the 
separation involves different regions of momenta. For very large bias the dispersion 
law reduces to the one-way hopping form exp(i ■ Pi/2), but multiplied by the sign of 
cos(pi) + cos(j>2)- this is the signature of the hopping disorder, which translates into a 
singular behavior along the Bragg lines. It is seen then that a perturbative approach 
fails also in the extremal regimes. 



5 Appendix Al 

We first point out here a relevant property of gaussian averages of time-ordered expo- 
nential operators, holding for perturbations delta-correlated in time. Let us consider 
the operator H(t) = Ho(t) + V(t), where the perturbation term is given through its 
correlator < V(t) (g> V(t') > = S(t - t') ■ A(t). One has, by definition: 

l=0,p=± Jti +1 =t' m=1 

The average leads to: 

<Texp[J dT-H(r)]> = Texp[ J dr ■ [Hq{t) + sign(t - t') ■ t4//(r)]] 
<V(t)-V(t')> = 2-S(t-t')-V eff (t) 



6 Appendix A2 

The calculation of the averaged transition probability (eq. |^), in the hermitean case, 
can be performed also by means of an exact resummation of ladder diagrams; the 
probability then coincides with the diffuson amplitude. The retarded and advanced 
Green's functions are: 

[-id t Ti r n + h{t)]-G ± {t,t') = 5{t-t') 

G ± (t,t') = ±i6{±(t-t'))-U(t,t') 

G+(t,t') = [G-(t',t)}+ 

< x\U(t,t')\x' > | 2 = < x\G+(t,t')\x' > | 2 + | < x\G-(t,t')\x' > | 2 



s 



The average of the U operator is performed first, obtaining an effective single particle 
generator h e ff(t): 



<U(t,t')> = Texp[- sign(t-t') f dr-ft e //(r)] 

Jt> 

heff{t) = ^2D(/j,;t)\x>< 



if 

Let us consider now the perturbative expansion for the average of < x\G + (t, t')\x' > 
■ < x'\G~ (t' ,t)\x >, where h(t) is the perturbation. One can show that the expansion 
can be written as the sum of particle-antiparticle diagrams where the 'free propagator' 
lines are substituted with the exact averaged Green's functions Gf v , and the contrac- 
tions are made only between particle and antiparticlc disorder vertices. In fact, when 
contracting particle-particle vertices one is computing contributions to the averaged 
propagator. Due to causality and to delta correlation in time, crossed diagrams are 
zero ; one is left with the sum of ladder diagrams, i.e. only the diffuson survives. The 
basic contribution to the ladder has the form: 

< t,x;t,x\M\t',x';t',x r >= J^(< x\h(t)\y > ■ < y\h(t)\x >) av 

y,y 

■ < y\G+ v (t, t')\x' > ■ < ~x*\G~ v (F , t)\y >= 5(t - t) ■ 5^ ■ S x , - ■ 6(t - t>) ■ 0(t - V) 

■ X! *) • l S x,x'- eii + <W+eJ ■ exp[(- ^2 ( / drD(jj,; r) - _ drD(fj,; t))] 

M ij, Jt ' •>*' 

Let us introduce the following operator: 

<t,x\P\t',x' >=6(t-t')J2 D (Kt)-[S x , X '-e ti +S x , x >+eJ-exp[-2Y^ f drD^r)], 

then use the identity: 

< t,x;t,x\M l \t',x';¥,x 7 >= S(t-t) ■ S x ^ ■ S x ,~ < t,x\P l \t',x' > . 
The diffuson amplitude A(i, x; t' , x') is given by: 

//"* °° 
dti6(t - ti) ■ exp[ - / dT ■ £>(/«; t)] < h^l^P^t^x' > 
n , ' tl 1=0 

= 6{t - t') < x\Texp[ [ dr ■ (S(t) -2V%;t))]|i / > 

S(t) = ^2D(fj,;t) ■ [\x >< x + e^ + \x + >< x\] 

x ,n 

In the momentum representation one recovers diffusion: 

A(t,x;t',x') = - / dk-exp[ik- (x-x') -2 V f dr ■ D{p,;r) ■ (I - cosk^)] 

(27T) J_ v J t i 
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7 Appendix A3 

Let us consider the interaction representation by taking the hopping term as a per- 
turbation; we have the evolution operator O: 

0(Mo) = ^2exp[- i [ dr ■ V(x; r)] \x >< x\ (19) 

The transformed hamiltonian h(t) is then: + (t, to) ■ ho(t) ■ O(t,to) = h(t); one 
verifies that h(t) is obtained from ho(t) through the substitution: 



u(x,[i;t) — > u(x f fi; t) • exp[i / dr ■ (V(x + e M ; t) — V(x; r))] (20) 

Jto 

w(x,fj,;t) —* w(x, fi;i) ■ exp[i / dr ■ (V(x + e M ; r) — V(x; r))] 



The invariance of correlators under this transformation is the origin of the indepen- 
dence of Eq. |^ on the potential. Notice further that the following identity holds in 
general: 

U Q (t,t') = d(t,t )-U{t,t')-d + (t',t Q ), (21) 

where Uq and U are the evolution operators of the original hamiltonian ho (see Eq. 
|l|) and of h respectively. 



8 Figure captions 

Fig. la: real part Er of E(P), over one half of the Brillouin zone (— tt < Pi < ir , 0. < 

P2 < tt), for the parameters D = 1., k = 0.5, a = 0.1 

Fig. lb: imaginary part Ei of E(P), same parameters as in Fig. la. 

Fig. 2a, b: competition between diffusion and bias: same as in Fig.l, but with D = 

1. ,A; = 0.6 ,a = 0.75 

Fig. 3a, b: strong bias:!) = 1. , k = 5. a = 0.75. 
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